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ABSTRACT 

The interacting binary white dwarf (AM CVn) systems HM Cnc and V407 have orbital 
periods of 5.4 min and 9.5 min, respectively. The two systems are characterized by 
an "on/off" behaviour in the X-ray light curve, and optical light curves that are 
nearly sinusoidal and which lead the X-ray light curves in phase by about 0.2 in both 
systems. Of the models that have been proposed to explain the observations, the one 
that seems to require the least fine tuning is the direct impact model of Marsh & 
Steeghs (2002). In this model, the white dwarf primary is large enough relative to the 
semi-major axis that the accretion stream impacts the surface of the primary white 
dwarf directly without forming an accretion disc. Marsh & Steeghs (2002) proposed 
that in this situation there could be a flow set up around the equator with a decreasing 
surface temperature the further one measured from the impact point. In this study, 
we estimate the light curves that might result from such a temperature distribution, 
and find them to be reasonable approximations to the observations. One unexpected 
result is that two distinct X-ray spots must exist to match the shape of the X-ray light 
curves. 

Key words: accretion, accretion discs — binaries: general, close — hydrodynamics — 
novae, cataclysmic variables — stars: individual(HM Cnc; V407 Vul) — white dwarfs 
— X-rays: binaries — X-rays: stars. 



1 INTRODUCTION 



The ultracompact binary stars V407 Vul ( RX J1914.4+2456, 
iMotch fe Haberllll995l:lMotch etal.ll 19961) and HM Cnc (RX 



J0806.3+1527. Ilsrael est al.lll999l : I Burwitz fe Reinschl feoOl) 
have been the subject of extensive study since their dis- 
covery. Believed to be members of the interacting binary 
white dwarf (AM CVn) class, each system displays only 
a single photometric period, presumed in both cases to 
be the orb ital period. For HM Cnc, P OI b = 321.5 s = 
5.36 min (|lsrael et al.1 Il999l : IStrohmaverl 120031). and for 
V407 Vul, Port = 569.4 s = 9.49 mi n (|Motch etafl 1 19961 : 
I Cropper et al. 19981: Strohmav er 2004). In a previous publi- 
cation ( Dolence. Wood, fe Silverll2008l . hereafter Paper I) we 
simulated the hydrodynamics of the accretion stream using 
the method of smoothed particle hydrodynamics (SPH). We 
found that the impact spot was small enough and thus the 
accretion stream ram pressure high enough for the specific 
kinetic energy of the flow to be advected and thermalized 
below the photosphere of the primary, such that it can be 
emitted "downstream" of the impact point in the equatorial 
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flow that should develop. In Paper I we discussed the phys- 
ical nature of that accretion-driven flow and suggested the 
general features of any plausible photospheric temperature 
distribution that should result. 

The observational constraints for the systems HM Cnc 
and V407 Vul were discussed in detail in Paper I (also see 
ICropper et al.ll2004T ). Briefly, the systems are characterized 
by "on/off" X-ray light curves that trail the opt ical light 
curves by roughly 0.2 in phase (|Barros et al.|[2007l . and ref- 
erences therein), no hard X-rays are detected, the observed 
periods are decreasing with time at a rate con sistent with 
gravi t ational radiation angular momentum loss l|Strohmaveil 
l2003l . 12004 iRamsav. Cropper, fe Hakalall2006l ). and flicker- 
ing is minimal or absent ( Barros et al.l 20071 ). 



There are several models that have been proposed 
to explain the observations of these two systems: inter- 
mediate polar, polar, unipolar inductor, and direct im- 
pact. First, in the intermediate polar model originally sug- 
gested for V407 Vul , the observed per i ods reflect a white 
dwarf spin period (IMotch et al.1 Il996l : Ilsrael et al.1 Il999l : 
iNorton. Haswell. fe Wvnnll2004 ). However, the lack of ob- 
served flickering, emission features, or any other observed 
(e.g., orbital) period, as well as the "on/off" character of 
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the X- ray light curve , argue against this model. Citing these 
issues ICropper et al.l (|l998l ) proposed that V407 Vul could 
be a helium analog of a polar (AM Her) system (see also 
iRamsav et aill2002bh . but the lack of observed He emission 
lines, polarization, and hard X-ray photons mak e it unlikely 
this model is correct for these systems (e.g., Ilsrael et al.l 
l2002h . 

The unipolar inductor (U I) model (|Wu et all l2002i ; 
iDaU'Osso. Israel, fc Stell3l2007h is the only model without 
a Roche- lobe filling secondary. Similar to the model pro- 
posed bv lGoldreich fc Lvnden-Belll (|l969h for the Io- Jupiter 
system in which the conducting moon Io generates an elec- 
tromotive force (EMF) as a result of the relative motion 
of the magnetic field of Jupiter. The weak plasma envi- 
ronment allows currents to flow, driven by this EMF, and 
results in observable bright spots where these currents are 
dissipated in the atmosphere of J upiter. The UI mode l was 
intriguing g i ven initial r eports bv lRamsav et all (|200Ct ) and 
Ilsrael et all |2003l . 120041 ) that the systems V407 Vul and HM 
Cnc (respectively) displayed X-ray and optical light curve s 
that were nearly anti-phased. However, iBarros et all (|2007l ) 
were able to show that in fact timing errors had crept into 
the previous optical/X-ray analyses, and that for both sys- 
tems the optical phase leads the X-ray phase by roughly 0.2. 
A roughly-constant phase offset suggests a common physical 
cause for both systems, but the corrected result presents a 
severe problem for the UI model. In the UI model it is ex- 
pected that the X-ray-bright footprints of the magnetic field 
lines that carry the currents should be nearly on the merid- 
ian that passes through the line of centres between the two 
starfl The expected UI geometry thus yields optical light 
curves anti-phased from the X-ray, assuming that the opti- 
cal modul ation comes f r om th e heated face of the secondary 
star. The IBarros et all (|2007j ) result that the optical leads 
the X-ray in phase by ~0.2 for both systems implies that 
the X-ray spots are rotated more than 90° away from the 
line of centres for both systems. First, it appears on the face 
of it that this cannot plausibly be a stable minimum en- 
ergy plasma current configuration in the co-rotating frame 
for both systems, but even if so, this geometry places the 
secondary star below the horizon of the X-ray spots thus 
removing the p resumed heat i ng so urce of the face of the 
secondary star. IBarros et all (|2005h discuss these and fur- 
ther difficulties with the model, and conclude it does not 
apply to these systems. 

The final propose d model is the direct impact model. 
iNelemans et all (|200ll ) showed that if initial Roche-lobe con- 
tact and the onset of mass transfer in the AM CVn systems 
occurs at a short enough orbital period, the accretion stream 
can impact the surface of the primary white dwarf directly, 
and thus not form a disc. iMarsh fc Steeghsl |2002f ) applied 
this to suggest that the observations for V407 Vul could be 
explained by just such an Algol-like model. Marsh & Steeghs 
(2002) proposed that the bulk of the accreted matter would 
thermalize below the photosphere, mix, and boil to the sur- 
face at a temperature cool enough to emit soft X-rays. The 
temperature excess should broaden in latitude as additional 
atmospheric material is entrained in the equatorial flow that 
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would be established, and the peak emission energy should 
shift toward the optical with distance from the upwelling 
point. 

Because no hard X-rays are detected from these sys- 
tems, for this model to be viable the accretion stream must 
have sufficient ram pressure to advect the specific kinetic 
energy of the accreted material below the photosphere be- 
fore th ermalizing. In Paper I w e simulated the system ge- 
ometry |MarsOTiIeeghi (|2002l ) proposed for V407 Vul and 
found using SPH simulations that this constraint is satis- 
fied. The density in the accretion stream is roughly Gaus- 
sian about the midpoint of the flow, with a 1-er half width 
of 23 km in physical units at the impact point. The con- 
tours of mass flux for the impact region on the surface 
of the primary white dwarf are well approximated by a 
normal bivariate (2-D Gaussian) distribution. Fitting the 
isoflux contours, we found = 164 km (in longitude along 
the equator) and ag = 23 km (in latitude). The 3-er el- 
lipse encloses 99% of the mass flux, and has an area of 
A = 7r(3 ■ 23)(3 ■ 164) km 2 = 1.1 x 10 5 km 2 or a frac- 
tion of the stellar surface of fa<j = 8.5 x 10~ J . The mean 
mass flux within this area is M/A ~ 600 g cm -2 s _1 for 
M = 10~ 8 M e yr" 1 , the mean density of the stream is 
p~2x 10~ 6 g cm -3 , and the velocity of the accreted mass 
is roughly v « y / 2GM 1 /R 1 ns 0.01c. 



2 THE MODEL 

Within the direct impact model, the location of the accretion 
stream impact spot is fixed in the binary co-rotating frame; 
the spot is on the accreting star's equator, and has an impact 
angle with respect to the line of centres that is a function 
of the mass ratio of the system and the radius (mass) of the 
primary white dwarf. For most physically-plausible system 
parameters yielding dynamically-stable mass transfer, the 
impact point has an impac t angle roughly in the range 80 
to 120° l|Barros et all 120071 ). In this geometry, the altitude 
of the accretion stream as viewed from the impact point is 
<45°, the bulk of the specific angular momentum is hor- 
izontal to the local stellar surface, and we expect that an 
equatorial flow will rapidly be established. 

The accreted material will penetrate far below the 
photosphere to a depth where the envelope pressure is 
of order the ram pressure. For the mass transfer rates 
typical of AM CVn syste ms of such short orb ital period 
(10 -8 < M < 10" 7 M Q vr- MBildsten et al.ll2006h , we expect 
fram ~ 10 13 dynes cm~ 2 , a pressure found at a depth 
of ~ 60 km ~ 0.006J ?! in a typic al O.6M0 white dwarf 
of T ea ~ 20,000 K jWoodl Il995h . In a non- interacting 
white dwarf of this mass and temperature, the tempera- 
ture, density, and degeneracy parameter at this depth are 
T ~ 600, 000 K, p ~ 0.1 g cm" 3 , and 77 « -4, respectively. 
The fractional mass of the layer exterior to this depth is 
Afsheii ~ 10 -9 M*, or roughly the mass accreted in a few 
days to a month at the M rates quoted above. If we assume 
instead that the stream penetrates to a depth of where the 
envelope pressure is an order of magnitude higher than the 
ram pressure, P cnv ~ 10 14 dynes cm -2 , we find the depth to 
be ~ 100 km or equivalently 10 -8 M* in mass. The tempera- 
ture, density, and degeneracy parameter at this depth in the 
unperturbed model are T ~ 1.2 x 10 B K, p ~ 0.6 g cm , 
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and rj « —3, r espectively. The observed ~65 eV X-ra y source 
temperature (jlsrael et alj 120031 ; IStrohmaverl |2008| ) is thus 
consistent with upwelling of material from these depths. 

The thermalized temperature of the accreted gas (~ 100 
keV) will be orders of magnitude higher - and the density 
orders of magnitude lower - than the surrounding environ- 
men t, leading to energetic co nvective overturn. In the orig- 
inal iMarsh fc Steegha (|2002h model, the authors proposed 
that the relative phasing of the optical and X-ray light curves 
might be explained by a temperature profile around the 
equatorial band characterized by a high temperature at the 
upwelling point and then a decreasing photospheric temper- 
ature as the angle from the impact point increases. Because 
the equatorial flow is in the same direction as the orbital 
motion, a distant observer will see the characteristic tem- 
perature increasing with time, culminating with the passing 
of the X-ray-bright spot across the face of the star - this is 
what causes the optical light curves to lead the X-ray light 
curves in phase. 

In the direct impact model the highest photospheric 
temperature should be some distance downstream of the ac- 
cretion stream impact point where the upwelling occurs, and 
the range in latitude over which the photospheric tempera- 
tures exceed the nominal mean stellar photospheric temper- 
atures should be of order but larger than the width of the 
accretion stream at impact. Using the results obtained in 
Paper I, the 3-cr profile enclosing 99% of the mass flux has 
a full width of ~138 km, or equivalantly 0.78° on a white 
dwarf of radius 0.0144 R e (M = 0.5 Af Q ). The length of the 
3-cr profile is 980 km, or equivalantly 5.6°. We expect that 
the further downstream from the impact point, the broader 
the zone should be in which the temperature is above the 
nominal mean photosphere, and the smaller that tempera- 
ture excess should be. 

We wrote a code in IDL that models a temperature dis- 
tribution on a cylindrical wall that represents the equato- 
rial band of the accreting white dwarf. We take the latitude 
limits of this cylinder to be ±15° for the calculations here. 
A schematic representation of the computational layout is 
given in Figure [1] The actual computational grid has 360 ele- 
ments in longitude (0° < (j> < 360° ) , and 100 elements in lat- 
itude (0° < 6 < 5°). Because of the symmetry with respect 
to 8 = 0°, we do not need to explicitly calculate the contri- 
butions for 9 < 0° since in our model T(<j>,-9) = T(<j>,0). 
Although integration over a spherical surface would yield 
slightly different results for very low inclination angles, it 
would not change the overall results that we present in this 
preliminary study. 

In our IDL code, we integrate numerically over the visi- 
ble grid elements (180 x 100) at each phase. We compute for 
each element the blackbody spectral energy distribution for 
that element's temperature, correct for projection angle (in- 
cluding inclination) and limb darkening, and then sum over 
all visible elements. Finally, the contribution from the rest 
of the stellar surface is included, and the integrated spectral 
energy distribution is convolved with the Sloan u'g'r'i' filter 
transmission functions for dir ect comparison with the ob- 
servational results compiled bv lBarros et all (|2007r ). For the 
X-ray band, we simply assume a window with 100% trans- 
mission from 20 to 60A - roughly the sensitivity band of the 
Chandra LETG spectrograph used bv lStrohmaverl l|2008h in 
his study of HM Cnc. 




Figure 1. A schematic view of the computational grid. The nu- 
merical integration is over the heavy grid represented. The grid 
with lighter lines is symmetric with respect to the equatorial 
plane, and does not need to be computed explicitly. 



Knowing that neither system is eclipsing, but otherwise 
uncertain about the system inclinations, we assume i = 45° 
in the model. At this inclination the cylindrical approxima- 
tion for the equatorial band is still valid. As the inclination 
decreases from 90°, the fraction of the visible surface made 
up by the equatorial band declines and the temperature con- 
trast or latitude limits of the band must increase to yield the 
observed light curve amplitudes. At very low inclinations the 
visibility drops for the hemisphere below the equator, and 
the lower latitudes of the equatorial band may be hidden (see 
Figures 2 and 3), again requiring an increased temperature 
contrast or increased latitude limits to recover the observed 
optical light curve amplitudes. 

There are countless models for the equatorial-band tem- 
perature distribution that could be employed. For example, 
we initially implemented a model in which there was a sim- 
ple linear temperature decline with distance from a compact 
high-temperature (X-ray bright) spot, but quickly discarded 
that model as it produced synthetic light curves with far 
too much X-ray luminosity over too great a range in longi- 
tude/phase. Our next model included a single compact high- 
temperature (X-ray emitting) spot. Downstream from this 
spot the temperature distribution is a decreasing function 
of both distance from the high-temperature spot, and from 
the equatorial symmetry plane. In this model the temper- 
ature falls off exponentially with longitude, and then falls 
off in latitude approximately as a Gaussian function that 
has a width that increases with longitude. We use the stan- 
dard kernel function of smoothed particle hydrodynamics to 
approximate the Gaussian, and let our temperature distri- 
bution within the latitude limits of the equatorial band be 
given by the following expression: 

T(4>,0) = T basc + (Ti - T^e-^^WiB, h e (</>)] (1) 
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where T baso is the assumed background temperature of the 
equatorial band of the primary star, and Ti is the peak tem- 
perature of the distribution (Ti ~ 65 eV). The temperature 
of the remainder of t he star is T*. The function W is the 
SPH kernel function (jMonaghan fc Lattanziolll985l '). which 
approximates a Gaussian distribution but is identically zero 
beyond twice the smoothing length h: 



) 2 + l 



w(e,h) = { |[2-(!)] ; 



0< r < 1, 



1 < f < 2, 



(2) 







> 2. 



Physically we expect that the thermalized energy from 
the accretion stream should spread out in latitude as a func- 
tion of distance from the upwelling point, and we assume the 
functional form 



ho 



< 



ho + (hi - h )4>/4>i (j>o < <i> < 0i 
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where in the results given below we use ho = 1.2° and hi = 
7.5°. 

Unlike our original linear-decline model, this one im- 
mediately gave results that reproduced the general observa- 
tional results, including the on/off X-ray pulse shape and 
the optical/X-ray light curve phase offset. However, as we 
discuss below, a model with a single high-temperature spot 
does not reproduce the width of the X-ray pulses, and there- 
fore we were forced to introduce a second X-ray spot ~65° 
upstream of the original spot, characterised by temperature 
To. So in practice we typically set (j>o = 65° and <f>\ = 260°, 
where the results are not terribly sensitive to the choice for 
01. 

Finally, because the e-folding length required to re- 
cover the observed phase offset is such that the temperature 
of a grid element obtained using Equation (1) extended be- 
yond 2-7T exceeds that for the same element directly, we set 
the temperature of each grid element to be the maximum of 
T{<j>,0) or T{4> + 2%,6). We show in Figures [2] and |3] contour 
plots of the temperature distributions used below to pro- 
duce synthetic lightcurves representing HM Cnc and V407 
Vul, respectively. 
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Figure 2. Contour plot of two cycles of the temperature distribu- 
tion representing HM Cnc. The base temperature is Tb aso = 100 
kK, and the contour levels are 100001 K, 100.5 kK, 110 kK, 150 
kK, and 200 to 700 kK in steps of 100 kK. The white lines plot- 
ted over the second cycle show the function hg(</>). The high- 
temperature feature at ~360° and ~720° is the impact point 
X-ray spot. 



o 

T3 
3 




200 400 
Longitude (deg) 



600 



Figure 3. Contour plot of two cycles of the temperature distribu- 
tion representing V407 Vul. The base temperature is Tb asc = 25 
kK, and the contour levels are 25010 K, 30 kK, 50 kK, and 100 to 
700 kK in steps of 100 kK. The white lines plotted over the second 
cycle show the function hg(<p). The high-temperature feature at 
~360° and ^720° is the impact point X-ray spot. 



3 THE OBSERVATIONS 



The p hase-folded light curves compiled by iBarros et al.l 
(|2007l ) are shown in Figure [3] The authors reanalysed the 
Chandra X-ray data for V407 Vul obtained on 2003 Febru- 
ary 19 (solid line) and 2003 November 24 (dashed line). 
The X-ray data f or HM Cnc were digitized and fit from 
IStrohmaverl (|2005T ). Barros et al. (2007) obtained the optical 
light curves using the ULTR ACAM triple-beam CCD cam- 
era (|Dhillon fc Marshll200lh on the 4.2-m William Hershel 
Telescope in 2003 May and on the VLT in 2005 November. 
In the cours e of their ana l ysis, t he authors discovered a likely 
error in the IStrohmaverl l|2004h V407 Vul eph emeris, and in 
the optical timings of Israel et al.l l|2003l . |2004| . see Barros et 
al. (2007) for full details). 



In both V407 Vul and HM Cnc, the X-ray light curves 
are characterised by a pulse shape that rises steeply from 
zero detected flux in the "off" state to maximum flux, and 
then falls more gradually from peak light. Additionally, in 
both systems the phase offset between the X-ray and opti- 
cal data is roughly 0.2 in phase, with the opti cal peak flux 
leadin g the X-ray. Contrasting the systems, the lBarros et al.l 
l|2007h data show that, whereas the fractional amplitudes of 
the V407 Vul optical light curves increase dramatically from 
the r' to v! bands, the HM Cnc fractional amplitudes are rel- 
atively constant across all optical wavelengths. In addition, 
whereas the V407 Vul optical light curves are very nearly 
pure sinusoids, the HM Cnc light curves reveal a second- 
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Figure 4. Phase-folded light curves of V407 Vul (left- hand panel) obtaine d 2003 February 19 (solid line) and 2003 November 24 (dashed 
line) and HM Cnc (right-hand panel) as presented in iBarros et al. | (1200711 and kindly sent to us by S. Barros. Except for the HM Cnc 
j'-band data, the error bars are comparable to the size of the points, and have been omitted. The error bars for the X-ray data are not 
available, but can be estimated from the scatter in the data. 



harmonic that has an amplitude 15-20% that of the first 
harmonic. The optical spectrum of V407 Vul is dominated 
by a late-G (or early-K) type stellar spectrum (ISteeghs et al.l 
bOQfj ). Barros et al. (2007) suggest this star is either a line- 
of-sight coincidence in a dense field, or that it is the distant 
member of a hierarchical triple system, but that it cannot 
be the direct cause of the optical or X-ray variations. What- 
ever the ultimate determination, the longer- wavelength light 
curves are increasingly diluted from the flux from this star, 
causing at least in part the decreasing amplitude that is ob- 
served in longer-wavelength filters. We note that we omit 
the i'-band light curves from the analysis below because of 
the increased noise in the case of HM Cnc, and because of 
the dilution by the G star in the case of V407 Vul. 



4 RESULTS 

The primary goal of this project was to explore whether 
a plausible surface temperature distribution consistent with 
the direct impact model could give rise to the light curves ob- 
served for HM Cnc and V407 Vul, and in particular, whether 
this model could yield the phase offset and pulse shapes ob- 
served. The results we present in this section are representa- 
tive of the model output, but we caution that the solutions 
are not unique. 

4.1 The X-Ray Light Curves 

We begin by exploring the geometry and nature of the source 
of the X-ray light curves in both systems. In our model, we 
first assumed a small X-ray-bright region presumably just 
downstream of the impact point, and that the temperature 
distribution fell off exponentially as a function of longitude, 
while broadening in latitude. In Figure [5] we show a sin- 
gle cycle of the the X-ray phase-folded light curves for both 



systems with a simulated light curve also shown with and 
without limb darkening included. For V407 Vul, we include 
only the 2003 February 19 data for clarity. The simulations 
shown here are of a single 65 eV (750 kK) hot spot on the 
grid described above. The dotted line was calculated assum- 
ing no limb darkening, the solid line was calculated assuming 
a linear limb darkening law I(n)/I(l) — 1 — u x (l — fj,), and 
the short-dashed line was computed with a quadratic limb- 
darkening law I(n)/I(l) — 1 — u x (l — A* 2 )) where fi = cos(6?) 
and where we assume u x = 1.0 in both relations. Even with 
this extreme limb darkening coefficient, it can be seen that 
for both systems the observed light curves rise more steeply 
than the synthetic curve. One possible effect not included 
that could give the observed steep slope is a stellar surface 
that is non-spherical. The hydrodynamic simula tions of di- 
rect im pact accretion between polytropic stars bv lMotl et al.l 
l|200l see, e.g., their Figure 1) show just such an effect. 

Figure [5] also demonstrates that the pulse width in both 
systems is wider than would be expected for a single X-ray- 
bright spot rotating into and out of view. Interestingly, both 
systems show a gentler slope during part of their decline, and 
that slope appears to be roughly the same for both systems. 
In HM Cnc, the final decline is observed to be as steep as the 
initial rise, which is to say as steep as a single bright spot on 
the stellar surface. To demonstrate this we copied a portion 
of the single-spot fully-limb-darkened curve and shifted it 
by 0.21 in phase, and have shown this in Figure [5] as the 
long-dashed line. We find that to match the general shape 
of the X-ray light curves, we were forced to include a second 
compact X-ray spot some 60 to 65° upstream of the primary 
X-ray spot and small enough that it does not contribute sig- 
nificantly to the optical flux. We speculate that the primary 
spot is the upwelling point, and the secondary spot is the 
region around the impact spot. While most of the specific 
kinetic energy of the accreted mass is advected and ther- 
malized below the photosphere, the low-density outer sheath 
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Figure 5. Phase- folded X-ray light curves (filled circles) of V407 
Vul (left-hand panel) and HM Cnc (right-hand panel). Also shown 
are synthetic X-ray light curves calculated with three different 
limb-darkening laws: no limb darkening (dotted line), a linear 
limb darkening law with a coefficient u x = 1.0 (solid line), and a 
quadratic limb darkening law again with u x = 1.0 (short-dashed 
line). Also shown is the light curve with the quadratic relation 
offset by 0.21 in phase (long-dashed line), demonstrating that the 
final decline in the X-ray light curve of HM Cnc is well matched 
by a second compact X-ray spot. 



of the accretion stream ( > 5-cr from the midline) will have 
Pram < f'photosphcrc and this should yield an X-ray emitting 
halo around the impact point. There will be a strong surface 
flow towards the impact spot in this region with speeds of 
order the sound speed as a result of the higher- density core 
of the accretion stream sweeping any and all atoms deep be- 
low the photosphere; effectively, the high-density core of the 
impact spot is a hole in the surface of the star into which 
the surrounding materi al will flow. The dir ect impact hydro- 
dynamic simulations of iMotl et al.l l|2007l ). while computed 
for a mass ratio of q = 0.4 and thus higher than thought to 
apply to HM Cnc and V407 Vul, do indicate an upwelling 
point substantially displaced from the impact point. 

4.2 HM Cnc 

Figure [5] shows a sample model result for HM Cnc plotted 
over the data. To determine suitable model parameters for 
each system we first find those that yield the best match "by 
eye" to the X-ray light curve, and then explore the remain- 
ing parameters to find a reasonable match to the optical 
light curves. The parameters for the HM Cnc model shown 
in Figure [5] are given in Table 1 and the contour plot of the 
temperature distribution is shown in Figure [2] but again we 
emphasize that this set of parameters is almost certainly not 
unique and should be taken as representative, not definitive. 
As an example of this, no set of parameters that we found 
could cleanly reproduce the excess flux in the optical be- 
tween phases 0.8 and 1.0, but by artificially doubling the 
temperature of each grid element between longitudes 220° 
and 260°, we obtain a result which is a better approxima- 
tion to the data. The physical source of this excess flux is not 
clear. If the broad outlines of the current model are correct, 
then the source cannot be the heated face of the secondary 
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Figure 6. HM Cnc data with a sample model result overplotted. 
The model reproduces the X-ray pulse shape, the phase offset, 
and the general character of the optical light curves. 



star, as the phasing is is almost exactly wrong. For an as- 
sumed impact angle of roughly 100° the upwelling point is 
at an angle of ~ 160 and the excess optical flux is at an angle 
of ~ 340°, or roughly when the back side of the secondary 
is towards t he ob server. The hydrodynamic simulations of 
IMotl et all (|2007l ) show some evidence for large amplitude 
standing waves on the surface of the accreting star after the 
outer layers are spun up and perhaps this could be related 
to the observed excess flux. Additional work will be required 
to reveal the source of this radiation. 

There are some physical conclusions we can make from 
the simulations that most closely approximate the data. The 
synthetic X-ray light curve is a reasonable approximation to 
the data, but the data rise more steeply to maximum, per- 
haps because the surface of the star is slightly non-spherical. 
Our experiments clearly show that two distinct X-ray spots 
are required to reproduce the observed X-ray pulse shape - 
a single extended high-temperature ridge results in slopes of 
rise and fall that are shallower than those observed. Obser- 
vationally, the optical light curves all have nearly the same 
fractional amplitude regardless of filter, implying that the 
underlying source must have a very high temperature that 
is fluctuating, such that the optical passbands are far out on 
the Rayleigh- J e ans ta il of the spectral energy distribution. 
iBildsten et all (2006) suggest that compressional heating 
should result in these stars having temperatures T e g > 100 
kK, but we find fractional amplitudes of only ~0.04 if we 
assume the entire stellar surface has a temperature of 100 
kK, and not just the equatorial band base temperature. 

4.3 V407 Vul 

Figure [7] shows a sample model result for V407 Vul, ob- 
tained using the parameters given in Table 1. The temper- 



Synthetic Direct Impact Light Curves 7 



Table 1. V407 Vul, HM Cnc Model Parameters. 



Parameter 


V407 Vul 


HM Cnc 


Description 


T 


750 kK 


880 kK 


Impact spot temp 


Ti 


750 kK 


750 kK 


Upwelling spot temp 


<£o 


60° 


65° 


Angle between spots 


fa. 


260° 


260° 


Temp distribution param. 


^"band 


25 kK 


100 kK 


Boundary layer base temp 


T* 


25 kK 


25 kK 


Stellar photosphere temp 


T 4> 


100° 


80° 


Temp e-folding angle 
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Figure 7. A sample model result plotted over the observations of 
V407 Vul. The model reproduces the amplitudes and shapes of the 
X-ray and u' light curves reasonably well (solid lines). The G-star 
spectrum dilutes the observed amplitudes in longer passbands. 



ature distribution is shown in Figure [3] The direct impact 
model provides a reasonable approximation to the obser- 
vations. As discussed above, the optical spectrum of V407 
Vul is dominated by a G-type stellar spectrum, so the am- 
plitude variations observed are not intrinsic to the direct 
impact scen ario, and as a r esult are not reproduced by the 
simulations. ISteeghs et al. (2006) note that the amplitudes 
are consistent with a constant G-star spectrum plus a vari- 
able blue component. For this model we show the g' and r' 
synthetic light curves as dotted lines. In order to recover the 
observed phase offset for this star, we must use an e-folding 
length of ~ 100°, but this then yields excess flux in the 
X-ray light curve in advance of the observed rise. 



5 CONCLUSIONS 

W e have shown tha t a var iant of the direct impact model 
of lMarsh fe Steeghsl (|2002l ) can reproduce the general fea- 
tures observed from the ultracompact binary AM CVn sys- 
tems HM Cnc and V407 Vul. These preliminary results are 
promising, but should be taken as representative only. Im- 



proved fits to the observations would be possible by iter- 
ating to a minimim \ 2 using differential corrections to the 
model parameters as is implemented in the Wil son-Devinney 
code I Wilson fc Devinnevl [l97ll ; I Wilson! Ii"994h . or to allow 
the temperatures of the current grid points to vary freely 
(with or without constr aints) and use a genetic algorithm or 
related technique (e.g.. lGoldberd fl989) to find solutions for 
the temperature distributions that match the light curves 
of the individual systems. The solutions found in this way 
would not be unique as the number of grid points exceeds 
the number of data points, but additional physical insight 
might be gained through suc h a study. This is similar to 
the "f irefly" approach us e d bv l Hakala. Cropper, fc Ramsavl 
|2002l ) and iBridge et al] (|2004F in their studies of eclipse 
mapping of polar accretion streams. 

The primary results of this work are 

(1) The direct impact model of lMarsh fc Steeghsl (|2002l ) 
is a plausible mechanism to produce the observed X-ray and 
optical light curves observed in the ultracompact AM CVn 
binaries HM Cnc and V407 Vul. Using a reasonable guess for 
the general parameters of a temperature distribution model, 
the optical/X-ray phase offset is reproduced, as well as the 
general shapes of the light curves themselves. Subtle differ- 
ences between the model results and the observations indi- 
cate that the model is still incomplete. 

(2) In both systems the X-ray light curves rise signifi- 
cantly more steeply than would be predicted for a single 
compact hot spot on the surface of a spherical star, even 
assuming an extreme limb-darkening law and coefficient for 
the X-ray band. We speculate that the surface of the star 
may be non-spherical as a result of the ongoing accretion, 
and that the upwelling X-ray region in particular may have 
a photospheric radius slightly larger than than the average 
for th e star. The hydrodynamic simulations of iMotl et al.l 
(2007) appear to support this scenario. 

(3) The X-ray light curve from system HM Cnc and V407 
Vul can only be reproduced by including a second X-ray- 
bright spot ~ 65° upstream of the primary (brighter) X-ray 
region. To fit the observed light curve, there must be two dis- 
tinct spots - a high-temperature ridge spanning this range in 
longitude will not yield an X-ray light curve that reproduces 
the observed pulse shape. We speculate that the secondary 
spot is the actual impact point of the accretion stream, and 
that the prim ary spot is the re gion where the upwelling oc- 
curs (again cf. lMotl et al.|[2007l ). The X-rays from the impact 
region most likely originate from the low density outer re- 
gions of the accretion stream for which P ra m ^ ^photosphere • 
We might expect that the X-ray emission from the turbu- 
lent upwelling point might be highly variable - considerably 
moreso than the X-ray emission from the impact point - and 
a detailed statistical analysis of the existing X-ray time se- 
ries might shed light on the viability of this proposed model. 

(4) For the accreti on rates predicted for these two systems 
(Bilds ten et al.l 2006), the ram pressure is sufficient to pen- 
etrate to a depth of order 100 km below the photosphere or 
10" 9 to 10" 8 A'i* in mass depth - i.e., log(l - M r /M*) « 8.5, 
where M r is the mass contained within radius r. The tem- 
peratures at these depths in the envelope of a 20 kK white 
dwarf model are characteristic of the 65 eV temperatures 
measured from X-ray spectroscopic observations. 

(5) The observed optical light curves of HM Cnc all have 
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roughly the same amplitude. To obtain this in the simu- 
lations, the base temperature of the equatorial band must 
be of order 100 kK, but the rest of the star must be at a 
lower temperature to re produce the observed fractional am- 
plitudes in the optical. iBildsten et al.l (|2006T > suggest that 
compressional heating of the accretor should lead to mean 
photospheric temperatures of T e g > 100 kK - marginally 
consistent with our simple model results. The amplitudes 
of the optical-band light curves from V407 Vul increase dra- 
matically from the i! band to the u' band. Be cause the opti- 
cal li ght is diluted with the G-star spectrum l|Steeghs et al.l 
2006) it is unknown how much of the amplitude trend in 
V407 Vul is the result of this dilution, and how much might 
be intrinsic to the accreting star's temperature distribution. 

(6) In order to recover the observed phase offset between 
the optical and X-ray light curves in V407 Vul, the model 
yields non-zero X-ray flux significantly in advance of the ob- 
served rise. This suggests that the actual temperature distri- 
bution resulting from direct impact accretion may be more 
complex than the simple functional form we have assumed 
for this preliminary study. 

While the current results are intriguing, what remains 
to be done is to model the hydrodynamics of the direct im- 
pact scenario directly as it applies to the envelope and atmo- 
sphere of the accreting white dwarf. Direct hydrodynamical 
simulations in 3D should allow the estimation of the photo- 
spheric temperature on the surface of the primary star, and 
from that it is straightforward to calculate an estimate of 
the phas ed light curves that s hould result. The hydrocode 
PLUTO dMignone et al.ll2007t ) is well suited to such prob- 
lems, and we are in the initial stages of implementing this 
project. We anticipate that explicit hydrodynamical simula- 
tions will not only reveal the true structure and dynamics 
of the outer layers of the accreting white dwarfs in HM Cnc 
and V407 Vu l, but also the nature of t he closely related sys - 
tem ES Cet dWarner fe Woudtll2002l ; lEspaillat et all 120051 ) 
where time-series spectroscopy reveals unusual features that 
may result from a truncated accretion disc, possibly result- 
ing from a grazing impact of the accretion stream with the 
white dwarf atmosphere as the system makes the transition 
from a direct impact system to a normal AM CVn accretion 
disc. 
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